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ABSTRACT 


Aft  Head  sewage  waste  generation  data 
collected  at  15-mlnute  Intervals  onboard  the  USS 
HAROLD  J.  ELLISON,  (DD  864)  both  in  port  and  at 
sea,  are  analyzed  by  Box-Jenklns  time  series 
analysis  techniques.  Results  show  that  both  in- 
port  and  at-sea  data  may  bt  modeled  by  pure  auto¬ 
regressive  Gaussian  Processes  with  dependence 
order  three  for  in-port  data  and  order  one  for 
at-sea  data.  Maximum  likelihood  parameter  esti¬ 
mates  are  provided  and  model  fit  checks  are 
performed.  A  method  is  provided  to  extrapolate 
model  results  to  other  U.S.  Navy  ships. 

ADMINISTRATIVE  INFORMATION 

This  task  was  accomplished  under  Work  Unit  1-2864-503  as  part  of  an 
overall  project  to  characterize  the  waste  streams  of  U.S.  Navy  ships  as 
tasked  to  this  Department  under  the  Naval  Environmental  Protection  Support 
Service  (NF.PSS)  NCBC  Project  Order  8-0011. 

INTRODUCTION 

The  David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(DTNSRDC)  surveyed  sewage  waste  generation  aboard  the  USS  HAROLD  J.  ELLISON 
as  part  of  an  ongoing  effort  to  characterize  the  non-oily  aqueous  waste 
streams  discharged  to  the  environment  by  U.S.  Navy  ships.  Shipboard 
waste  surveys,  including  this  particular  effort,  are  generally  sponsored 
by  the  Navy  Environmental  Protection  Support  Service  (NEPSS).  In  addition 
to  the  usual  quantity  and  quality  monitoring  performed  over  the  period 
5  November  1974  to  23  March  1975,  specific  waste  sources  were  monitored 
continuously  over  shorter  periods  to  provide  Information  on  flow  rate 
fluctuations  at  small  time  Increments. 

This  report  describes  the  methodology  used  in  analyzing  the  flow  rate 
fluctuation  data  obtained  from  continuous  monitoring  of  the  Aft  Enlisted 
Head  (combined  conmode  and  urinal  flush  wastes)  and  the  resulting  random 
process  models  deduced  from  the  analysis.  Although  other  major  waste 
sources  onboard  the  ELLISON  were  also  monitored  continuously,  this  report 
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is  specifically  limited  to  the  Aft  Enlisted  Head  source  to  demonstrate 
methodology  and  approach. 

In  the  following  report,  it  is  assumed  that  the  Aft  Enlisted  Head  flow 
rate  fluctuations  can  be  characterized  by  a  specific  form  of  the  family  of 

Autoregressive  Integrated  Moving  Average  (ARIMA)  models  described  by  Bo* 

1* 

and  Jenkins.  The  basic  philosophy  and  techniques  of  the  Box-Jenkins 
stochastic  model  building  approach  (Identification,  estimation,  checking) 
are  used  in  this  report,  supplemented  by  more  precise  statistical 
techniques  in  Che  model  identification  portion  of  the  analysis. 

BASELINE  FLOW  DATA  DESCRIPTION 

The  baseline  flow  data  to  which  the  analysis  will  be  applied  were 
obtained  during  the  period  2-12  December  1974  from  the  Aft  Enlisted  Head 
soil  drains.  The  recording  method  was  an  event  recorder  which  provided  an 
Indicator  mark  on  a  strip  chart  each  time  a  total  of  10  gallons  of  waste 
passed  through  a  metering  device.  Indicator  marks  wore  counted  in  each 
15-minute  period  during  the  ten-day  period.  Small  gaps  (usually  les9  than 
one  hour)  occurred  at  the  time  of  strip  chart  replacement  but,  due  to  their 
Infrequent  occurrence  and  short  duration,  were  assumed  to  represent  zero 
flow  periods.  The  resulting  data  provided  a  continuous  flow  record  for  a 
three-day  at-sca  period  (3-5  December)  and  a  six-day  ln-port  period  (6-11 
December).  Preliminary  evaluation  of  the  data  Indicated  that  6,  7,  and  8 
December  had  significantly  different  flow  characteristics  from  the  rest  of 
the  period  because  of  the  inclusion  of  a  debarkation  day  (6  December)  and 
a  weekend.  These  data  were  subsequently  eliminated  from  consideration  in 
the  analysis,  leaving  three  full  days  of  waste  flow  data  for  normal  ln-port 
and  at-sca  operations. 

The  data  presented  in  the  Appendix  are  in  units  of  tens  of  gallons 
per  15  minutes.  Figures  1  and  2  show  in-port  and  at-sea  flow  fluctuations, 
respectively,  for  one  24-hour  period. 


*A  complete  listing  of  references  is  given  on  page  25. 
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Figure  1  -  Sample  Flow  Fluctuation 
IN  PORT 
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10  GALLONS  PER  15  MIN 


Z  O 


MODEL  DEVELOPMENT 


GENERAL  MODEL  DESCRIPTION 

The  family  of  ARIMA  random  process  models  described  by  Box  and 
Jenkins  provides  a  mathematical  representation  of  a  wide  variety  of 
characteristics  of  random  fluctuations  observed  in  time  series.  Stationary 
(constant  mean),  linear  non-stat ionary  (random  mean  level  changes),  and 
seasonal  (cyclic-pattern  repetition)  are  easily  represented  in  this 
family  of  models. 

Since  the  subsequent  models  identified  are  of  the  stationary  form,  we 
will  briefly  summarize  the  characteristics  of  this  type  of  model.  If 
z^,  t“l,2,...  represents  the  random  values  associated  with  a  stationary 
time  series,  then  the  model  form  is 


z 

t 


c  *  Vt-l  *  Vt-2  * 


+  $  z 
P  t-p 


+  at  -  Vt-l  •  Vt-2 


6  a 
q  t-q 


where  c  is  a  constant 


are  constant  coefficients 

a  is  a  Normal  (Gaussian)  random  variable  with  zero  mean  and 


variance  o 

a 

p,q  are  the  dependence  orders  associated  with  and  , 
respectively 

As  the  model  indicates,  the  current  value  of  z ^  is  dependent  on  propor¬ 
tional  values  of  p  previous  z's  and  q  previous  a ' s .  The  dependency 
relations  are  not  intended  to  display  causal  relationships  but  merely 
reflect  the  fluctuation  pattern  as  z  evolves  in  time.  No  physical  inter¬ 
pretation  can  therefore  be  assigned  uoless  additional  information  is 
available. 

The  model  Is  developed  by  first  identifying  the  values  p  and  q  which 
define  the  order  dependence,  estimating  the  parameter  values  t ^ , 
1-1,2,. ...p,  8^.  J-l , 2 , . . . , q ,  and  then  checking  the  fit  by  analysis  of 
residual  values.  This  sequence  of  procedures  revolves  around  the  auto¬ 
correlation  estimates  and  defines  the  basic  Box- Jenkins  model  building 
philosophy.  The  analysis  procedures,  the  results  of  which  are  described 


below,  have  been  applied  to  the  "at-sea"  and  "in-port"  portions  of  the 
data  of  the  Appendix  separately. 

IDENTIFICATION  OF  MODEL  FORM 


The  autocovariance,  autocorrelation,  and  spectral  density  estimates 
for  the  first  twenty  lags  are  given  in  Tables  1  and  2  for  in-port  and 
at-sea  data,  respectively.  Figures  3  and  4  provide  the  corresponding 
information  graphically.  The  structure  of  the  flow  variation  shown  in 
Figures  1  and  2  and  the  pattern  of  the  autocorrelation  and  spectral 
density  estimates  Indicate  that  both  the  in-port  and  at-sea  fluctuations 
are  stationary  about  a  constant  mean  and  are  most  likely  pure  auto¬ 
regressive  (AR)  processes  (i.e.,  p  >  0,  q-0  in  the  general  model  described 
in  the  preceding  section  of  this  report). 

To  verify  these  initial  indications  and  obtain  optimal  values  for  p, 

2 

both  sets  of  data  were  passed  through  IMSL  subroutine  FTCOMP .  FTCOMP  did 
Indeed  verify  the  initial  estimates  and  provided  optimal  values  of  p-1  for 
the  at-sca  fluctuation  and  p**3  for  the  in-port  version,  with  q»0  for  both 
sets.  In  addition,  FTCOMP  confirmed  the  statlonarlty  assumptions 
previously  made.  These  results  imply  that  the  data  are  best  fitted  by 
pure  autoregressive  (AR)  models  of  the  form 
in  port 


z 

t 


c  *  *1*1-1  *  *2*t-2  *  * J*t-3  *  \ 


at  sea 


c  +  Vt-1  +  at 


These  models  may  be  designated  at  AR(3)  and  AR(1),  respectively. 
ESTIMATION  OF  PARAMETERS 

In  addition  to  optimal  selection  of  the  model  form,  FTCOMP  also 
provides  Maximum  Likelihood  Estimates  (MLE's)  of  the  parameters  of  the 
selected  model.  Table  3  shows  the  results  of  these  computations  for  both 
in-poi L  and  at-sea  models  and  also  shows  the  mean  and  variance  estimates 
of  each  time  series. 
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Figure  }  -  Autocorrelat ion  and  Spectral  Density  Estimates 

IN  PORT 


Figure  4  -  Autocorrelation  and  Spectral  Density  Estimates 

AT  SEA 


TABLE  3  -  MAXIMUM  LIKELIHOOD  ESTIMATES  OF  MODEL  PARAMETERS 


Mean* 

Var.** 

C* 

♦t 

o  2 

a 

At  Sea 

3.54 

23.72 

2.649 

♦j  -  0.252 

22.21 

In  Port 

1.32 

2.33 

0.596 

-  0.278 

♦2  -  0.109 

-  0.161 

1.95 

•In  units  of  10  *  gallons. 

••In  units  of  (10_1  gallons)2. 

MODEL  IDIOSYNCRASY 

ARIMA  modeling  requires  the  assumption  of  a  Gaussian  distribution  for 
the  white  noise  term  a£  and,  from  theoretical  considerations,  the  multi¬ 
variate  Gaussian  distribution  of  the  z's.  As  a  consequence,  under  certain 
conditions  on  the  mean  and  variance  of  the  distribution,  large  negative 
values  have  large  probability  of  occurrence.  This  situation  exists  for 
both  models  (in  port  and  at  sea)  and  is  incompatible  with  the  physical 
reality  of  non-negative  flow  rates.  To  compensate  for  this  idiosyncrasy 
in  simulation  applications,  the  following  function  should  be  used  with  the 
mod  els: 

i  z  if  z  -  0 
)  t  t 

w  -  < 

1  1 0  otherwise 

In  probabilistic  analyses,  the  usual  Gaussian  density  integrals  should  be 
adjusted  by  multiplying  the  results  by  l/(l-Prob(Zt  <  0))  where 

!0.19  in  port 

0.23  at  sea 
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DIAGNOSTIC  CHECKS  ON  FIT 

Two  types  of  diagnostic  checks  were  performed  to  determine  how  well 
the  model  fits  the  data.  The  first  type  checks  the  autocorrelation  of 
residuals  for  significant  residual  dependencies  and  the  second  checks  the 
spectral  density  for  suspicious  periodicities. 

Figures  5  and  6  are  graphic  displays  of  autocorrelation  and  spectral 
density  estimates  of  the  residuals  shown  In  Tables  4  and  5.  Under  the 
assumption  that  the  residuals  are  white  Gaussian  noise  with  variance  o^, 
then  the  autocorrelation  estimates,  r^,  of  the  residuals  at  lag  k  should 
not  exceed  ?2o(r^)  i  t2//n  (n  Is  the  number  of  sample  points)  for  k>0. 

Since  n  -  288  for  both  In  port  and  at  sea,  i2o(r^)  1  '0.118.  It  can  be 
seen  in  Tables  4  and  5  that  none  of  the  autocorrelations  of  the  residuals 
exceed  these  bounds,  except,  of  course,  at  lag  0.  A  more  formal  check  can 
be  performed  by  computing  the  statistic 

K  2 

Q  ‘  " 

k-1 

which  is  Chi-Square  distributed  with  (K-p)  degrees  of  freedom  (d.f.),  where 
K  is  the  number  of  lags  for  which  the  autocorrelation  estimates  have  been 
computed.  For  the  20  lags  shown  in  Tables  4  and  5 

Q  (in  port)  -  0.0)1  with  (20-3)  -  17  d.f. 

Q  (at  sea)  -  0.012  with  (20-1)  -  19  d.f. 

both  of  which  pass  the  0.25  tabulated  Chi-Square  values  by  substantial 
margins.  On  the  basis  of  these  diagnostic  tests,  both  models  provide 
excellent  fit  to  the  data. 

In  the  second  type  of  check,  the  spectral  estimates  of  the  residuals, 
the  spectral  densities  displayed  In  Figures  5  and  6  indicate  some 
possibility  of  additional  periodic  components.  However,  these  periodic 
components  are  thought  to  be  the  result  of  the  truncation  effects  described 
in  the  previous  section.  To  substantiate  this  assumption,  simulation 
results  for  each  truncated  model,  each  simulation  of  length  288,  /ere 
passed  through  the  spectral  analysis  ano  are  displayed  In  Tables  b  and  7 
and  Figures  7  and  8.  These  spectral  estimates  compare  favorably  In  Tables 
1  and  2  and  Figures  3  and  4  of  the  original  data.  Consequently,  It  mav  be 
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Figure  fe  -  Autocorrelation  and  Spectral  Density  Estimates  of  Residuals 
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Figure  7  -  Autocorrelation  and  Spectral  Density  Estimates 
of  In  Port  Simulation 


Figure  S  -  Autocorrelation  and  Spectral  Density  Estimates 
of  At  Sea  Simulation 


concluded  that  the  truncation  function  produces  the  apparent  periodicities 
indicated  in  the  spectral  estimates  of  the  residuals  and  no  expansion  of 
model  forms  is  required. 

MODEL  EXTRAPOLATION  TO  FLEET 

Although  the  flow  fluctuation  models  developed  In  the  previous 
sections  are  sufficient  to  characterize  the  Aft  Head  sewage  flow  rate 
onboard  the  USS  HAROLD  J.  ELLISON,  these  models  would  be  of  little  use 
unless  some  method  is  available  which  would  allow  application  to  other 
ships  of  the  U.S.  Navy  Fleet.  This  section  provides  such  a  methodology 
and  states  the  assumptions  under  which  this  extrapolation  Is  valid. 

The  Aft  Head  source  data  on  which  the  models  are  based  consist  of  the 
volumes  of  combined  commode  and  urinals  flush  water.  Since  flush  water 
is  directly  related  to  complement,  shipboard  crew  size  information  should 
be  the  principal  extrapolation  factor,  assuming  that  crew  habits  and  the 
quantity  of  water  released  per  flush  are  consistent  on  a  fleetwide  basis. 

A  standard  practice  is  to  repai r/replace  poorly  maintained  flushometers 
prior  to  the  start  of  a  NEPSS  Shipboard  Survey  so  that  a  consistent  flush 
volume  can  be  maintained  and  valid  comparisons  can  be  made  among  surveyed 
ships. 

The  second  assumption,  the  fleetwide  consistency  of  the  crew's  daily 
sanitary  habits,  is  presumed  to  be  reflected  in  the  autocorrelation 
function  pattern  of  the  data.  On  a  fleetwide  basis  in-port  sewage 
generation  is  represented  by  an  AR(3)  process  with  the  1*1,2, 3  values 

as  displayed  In  Table  3,  and  at-sea  generation  by  #\R(1)  with  *  0.252. 
Crew  habits,  however,  are  difficult  to  verify  and  consistency  must  be 
presumed  for  extrapolation  purposes. 

Under  these  assumptions,  the  only  adjustments  required  to  the 

,  .  2 

parameter  values  of  the  models  are  for  c  and  o^.  These  two  parameters  can 

he  made  dependent  on  ship's  population  level  In  the  following  manner:  c  and 

2 

0“  are  related  to  MEAN  and  VAR  of  Table  3  by  the  formula 
a 

c  -  MEANx(l 

12  p 

o2  -  VARx( 1  -  ♦  r  -  ♦,r,  -•••-♦  r  ) 
a  112  2  p  p 
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where  the  r^  are  functions  of  the  *j'«*  only.  MEAN  and  VAR  values  obtained 
from  the  Aft  Head  of  the  ELLISON  may  be  reduced  to  15-mlnute,  per  capita 
values  by  the  method  shown  In  Table  8.  The  results  of  Table  8  rest  on  the 
assumption  that  the  ratio  of  Aft  Head  flow  to  Total  Salt  Water  flow  Is 
equivalent  to  the  ratio  of  Aft  Head  users  to  Total  Ship  population.  This 
assumption  is  plausible  since  the  total  salt  water  flow  consists  of  ail 
commode  and  urinal  flow  on  the  ELLISON  (including  Forward  Crew's  Head  and 
Officer's  Head)  as  measured  by  survey  over  a  5-month  period.  It  Is  also 
assumed  that  there  is  Inconsequential  head  area  crossover  usage  between 
Aft  Crew,  Forward  Crew,  and  Officer  Quarters.  As  a  consequence,  under  the 
assumptions  stated,  designated  model  parameters  can  be  adjusted  to  reflect 
changes  In  ships'  force  level. 

To  demonstrate  the  procedure,  model  parameters  will  be  extrapolated  to 
reflect  the  sewage  generated  by  total  crew  of  the  ELLISON. 

With  an  Average  Total  Ship’s  Crew  of  135  In  port  and  205  at  sea,  we 

have 
At  Sea 

MEAN  •  205  x  0.0268  -  5.494  x  10  1  gallons/15  min. 

VAR  -  (205) 2  x  13.6823  x  10‘4  -  57.500  x  ( 1 0~ 1  gallons/15  min.)2 
C  -  MEANx(l-«  )  -  5.494(1-0.252)  -  4.110  x  10-1  Rallons/15  min. 

a2  -  VARx(l-*2)  -  57.500  x  l 1- < . 252) 2 ]  -  53.849  x  (10"1  gallons/ 

*  1  5 

15  min.)*1 


In  Port 

MEAN  -  135  x  0.01517  -  2.048  x  10'1  gallons/15  min. 

VAR  -  (135)2  x  3.1708  x  10'4  -  5.779  x  (10-1  gallons/15  min.)2 
C  •  MEANxd-^-H,-^)  -  2.048  x  (1-0.278-0.109-0.161)  - 

0.930  x  10“1  gallons/15  min. 

o:  -  VARx(l->,  r  -♦,r,-$,r,J)  -  5. 779 1 1- (0. 2  78)  (0. 349)-(0. 109)  (0. 206) 
a  112  2  3  3 

-(0. 161) (.256) ]  -  4.850  x  (10_1  gallons/15  min.)2 


where  Tj  • 

r2  •  Vl  +  *2 


»lr2  ♦  Vi  +  *3 
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TABLE  8  -  PER  CAPITA  MEAN  AND  VAR  ESTIMATES 
FOR  MODEL  APPLICATION 


Aft  Head  Gallons/Day 
Total  Saltwater  Gallons/Day 
Ratio 

Average  Total  Ship's  Crew 

2 

Estimated  Aft  Head  Users 
15  minute  MEAN  (TABLE  3) 

15  minute  VAR  (TABLE  3) 

15  minute  MEAN/CAPITA3 
15  minute  VAR/CAPITA" 


Averages  from  5  month  survey. 

Ratio  *  Average  Total  Ship's  Crew. 

15  minute  MEAN/Est imated  Aft  Head  Users 
(In  units  of  10  3  Ga 1 lons/Capi ta) 

2 

15  minute  VAR/ (Est Imated  Head  Users) 

In  units  of  (10  3  Gallons/capita)" 


In  Port 

At  Sea 

1664.9 

3181.0 

2587.8 

4929.4 

0.64  3 

0.645 

135 

205 

87 

132 

1.  32 

3.54 

2.40 

23.84 

0.01517 

0.0268 

3.1708  *  10_4 

j 

13.6823  *  10~4 
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